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We study the inverse energy transfer in forced two-dimensional (2D) Navier-Stokes 
turbulence in a doubly periodic domain. It is shown that an inverse energy cascade that 
carries a nonzero fraction of the injected energy to the large scales via a power-law energy 
spectrum cx k~ a requires that a > 5/3. This result is consistent with the classical theory 
of 2D turbulence that predicts a fc -5 / 3 inverse-cascading range, thus providing for the 
first time a rigorous basis for this important feature of the theory. We derive bounds 
for the Kolmogorov constant C in the classical energy spectrum Eik) — Ce 2 / 3 fc -5 / 3 , 
where e is the energy injection rate. Issues related to Kraichnan's conjecture of energy 
condensation and to power-law spectra as the quasi-steady dynamics become steady are 
discussed. 



1. Introduction 

It is well known that the advective nonlinearities of two-dimensional (2D) Navier- 
Stokes (NS) turbulence predominantly transfer energy to low wavenumbers (large scales) 
and enstrophy to high wavenumbers (small scales). The extent of this preferential transfer 



has been a subject of intense research since Fj0rtoft (1953) first noticed this interesting 



property of 2D turbulence. In the late 1960s Kraichnan (1967) postulated that for un 



bounded fluids in the limit of high Reynolds number, this preferential transfer achieves 
an extreme limit, by transferring virtually all energy to ever-lower wavenumbers (inverse 
energy cascade) and virtually all enstrophy to a high wavenumber k v 3> s (direct enstro- 
phy cascade) if the turbulence is driven by sources at intermediate wavenumbers around 
s. The transfer of energy and enstrophy in this extreme manner is known as the dual 
cascade. The energy is predicted to cascade via a fc -5 / 3 energy inertial range, and the 
enstrophy is predicted to cascade via a fc~ 3 enstrophy inertial range. The dissipation 
wavenumber k v determines the region where the enstrophy gets dissipated. The resulting 
dynamics is quasi-steady since the inverse energy cascade presumably proceeds indefi- 
nitely in time toward wavenumber zero. This theory was later advanced by Leith (1968)[ 



Batchelor (1969) and Kraichnan (1971) and has become a classical theory of 2D turbu- 



lence, known as the KLB theory. 

For a fluid confined to a doubly periodic domain, the inverse energy cascade is halted 
at the lowest wavenumber fco corresponding to the integral length scale of the system. 
This arrest of the inverse cascade may be met with complex responsive adjustments 
of the fc -5 / 3 range. In any case a rise of the total energy in the available energy range 
necessarily occurs. According to Kraichnan (1967)) such a rise occurs only at fco, resulting 



in what can be termed an 'energy condensate' that singlehandedly carries most of the 
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total kinetic energy. f Assuming that the energy spectrum eventually becomes statistically 
steady, with a persistent inverse cascade up to ko , the energy condensate not only would 
carry virtually all the system energy but also would account for virtually all the energy 
dissipation. This means that Kraichnan's energy condensate would also be an enstrophy 
condensate, although the latter would be of a lesser degree. For a sufficiently wide energy 
inertial range, the enstrophy dissipation by the energy condensate would be negligible. 
This possibility allows for the enstrophy dynamics to be virtually unaffected by the 
condensate and for the direct enstrophy cascade to remain intact in Kraichnan's picture. 

Numerical results concerning the Kraichnan condensate are inconclusive and, in fact, 
controversial. Smith & Yakhot (1993, 1994) argue for a positive answer to the Kraich- 
nan conjecture. On the contrary Borue (1994) and Tran & Bowman (2004) observe that 
after the arrest of the inverse energy cascade, growth of energy and enstrophy occurs 
throughout the energy inertial range. More precisely they find that as the turbulence 
approaches a steady state, a k~ 3 energy spectrum forms at large scales. Although this 
spectrum still means that most of the energy is concentrated at the lowest wavenumbers, 
it is by no measure close to the Kraichnan picture. In particular a k~ 3 energy spectrum 
means a k~ x enstrophy spectrum, so that no accumulation of enstrophy occurs at the 
largest scales and the dissipation of energy cannot occur primarily at or in the vicinity 
of ho- 

Theoretical studies have shown under a variety of assumptions that for a doubly peri- 
odic domain, the ratio of enstrophy dissipation to energy dissipation equals s 2 — the ratio 
of enstrophy injection to energy injection (Constantin et al. 1994; Eyink 1996, p. 110; 
Tran & Shepherd 2002; Kuksin 2004). These authors have concluded from this result 
that there can be no direct enstrophy cascade. However this inference assumes that all 
the enstrophy is actively involved in the enstrophy cascade, an assumption that is in- 
consistent with the Kraichnan conjecture. Indeed, several of these studies have explicitly 
or implicitly excluded a priori the Kraichnan conjecture in favour of an assumption of 
power-law spectra. There is thus a need to explicitly examine the extent to which the 
Kraichnan conjecture is consistent with analytical constraints on 2D turbulence. 

In this paper we study several aspects of the KLB theory. In § 2 we introduce some 
necessary notation and derive a basic identity and inequality. In § 3 we give a general 
review of the KLB theory, including a discussion of the dual-cascade hypothesis and of 
the Kraichnan conjecture of energy condensation. In § 4 we derive an upper bound for 
the inverse energy flux and show that a nonzero steady flux requires that the inverse- 
cascading range be at least as steep as k~ 5 ^ 3 . More precisely, we show that an energy- 
transfer range shallower than fc~ 5 / 3 is unable to support a constant inverse energy flux. 
This result provides for the first time a rigorous basis for the Kolmogorov-Kraichnan 
fc -5 / 3 inertial range. We also derive bounds for the Kolmogorov constant C in the classical 
energy spectrum E{k) — Ce 2 / 3 fc -5 / 3 , where e is the energy steady injection rate. In § 5 
we derive a constraint on the enstrophy and palinstrophy distribution in the Kraichnan 
picture and contrast it with power-law scalings. Concluding remarks and discussion are 
given in § 6. 



f Kraichnan (1967)] also predicts that as the energy accumulates at ko, the —5/3 spectrum 
is modified toward absolute equilibrium of the form E(k) tx k/(/3k 2 + a), where j3 and a are 
constant (see Kraichnan 1967, p. 1423a). Note that Kraichnan's conjecture would be trivially 
extended to an energy condensate in the vicinity of ko, not just at ko- In this work we address 
this more general scenario. 
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2. Mathematical preliminaries 

In the vorticity formulation the forced 2D NS equation governing the motion of an 
incompressible fluid confined to a doubly periodic domain [0,L] x [0,L] is 

d t AV> + JO, AVO = v& 2 ip + /, (2.1) 

where ^(x, t) is the streamfunction, J ~(6, *&) — 9 x d y — 9 y i) x , v the kinematic viscosity and 
f(x,t) the forcing. The velocity field v(x,t) can be recovered from the streamfunction 
ip(x,t) by v = (—tjjy,ip x ). The nonlinear term admits the identities 

(MM)) = -<M&0)) = -<MM)), (2-2) 

where ( • ) denotes a spatial average. As a consequence we have the twin constraints 

(MV>, A ^)> - = (AMV>, AV)), (2.3) 

so that the energy £7 = (|V?/'| 2 )/2 and enstrophy Z = (|A-0| 2 )/2 are conserved by 
nonlinear transfer. 

We now derive a simple identity and an inequality, which are used in § 4 to calculate the 
nonlinear triple-product term representing the inverse energy transfer. By straightforward 
calculation we have 

A J[6, i?) = J(A9, z9) + J(9, Atf) + 2J{6 x ,d x ) + 2J{9 y , <& y ). (2.4) 

Hence 

(A0J(O, 0)) = (0A J(6, ??)) = (0J(A0, i?)) + 2(9J(9 X ,# X )) + 2(9J(e yi 4 y )) 

= -(AflJ(M)) + 2(9J(9 X ,# X )) +2(0J(0 y ,dy)). (2.5) 

It follows that 

(A0J(M)) = (9J{9 x ,ti x )) + (M^ir))- (2-6) 

For the present case of a doubly periodic domain of size LxL, the Fourier represen- 
tation of ip{ x ) is 

V>(x) = ex P{* fc ' x}^(fe). (2.7) 

k 

Here fe = fco(ji, to), where fco = 2n/ L is the lowest wavenumber and n and to are integers 
not simultaneously zero. For a given wavenumber I let - < and %j> > denote, respectively, 
the components of if) spectrally supported by the disk d = {k : k < £} and its complement 
D = {k : k > £}, i.e. 

'0 < = cxpjzfc • x}^(fe), ifr* = /J expjifc • x}tp(k). (2-8) 

feed fcG-D 

The lower-wavenumber component tp < satisfies 

1/2 , v 1/2 

fc 4 |^(fe)| 2 = C ^< /2 , (2.9) 



|A^<|<^fc 2 |^(fc)|< (fel) 

fcGd \fcGd / \fcerf 



where c is an absolute constant of order unity and Z < = (|AV' < | 2 )/2 the large-scale 
enstrophy density associated with the wavenumbers k < i. In (|2.9() the Cauchy-Schwarz 
inequality is used in the second step, and the sum J2ked 1 ~ ^ 2 Ao represents the number 
of wavevectors in d. 



4 



C. V. Tran and T. G. Shepherd 



3. The KLB theory 

This section reviews the central features of the KLB theory: the dual-cascade hy- 
pothesis and the Kraichnan conjecture of energy condensation. The discussion deviates 
significantly from the original works that lead to the dual-cascade hypothesis, thereby 
providing a new look at the KLB theory. We also discuss some recent results, which, on 
the one hand, can modify KLB significantly and, on the other hand, cast doubt on some 
aspects of the theory. 



We consider turbulence driven by sources localized around a wavenumber s that supply 
an energy injection e > and an enstrophy injection s 2 e. Let e(fe) and k 2 e(k) be nonlinear 
redistributions of these injections in wavenumber space. Given arbitrary wavenumbers 
r\ < s and r 2 > s we have 



where the second inequalities in l|3.1|) and (|3.2I) are, respectively, due to the conservation 
of enstrophy and energy. It follows that the respective fractions of e and of s 2 e that 
get transferred to k > r 2 > s and to k < r\ < s are bounded from above by s 2 jr\ 
and by r 2 /s 2 , respectively. These constraints correspond to the well-known prohibition 
of a significant direct (inverse) transfer of energy (enstrophy) : no considerable fraction 
of e (s 2 e) is allowed to get transferred to wavenumbers k ^> s (k <C s). Hence when 
these injections spread out in wavenumber space, most of the energy (enstrophy) gets 
transferred toward lower (higher) wavenumbers, a qualitative conclusion that is behind 
the dual-cascade hypothesis. This dynamical behaviour of energy and enstrophy is well 
confirmed by numerical simulations reported in the literature. Note that the preferential 
transfer described above reflects the collective effects of all admissible nonlincarly inter- 
acting triads. Detailed analyses of individual triads may not lead to the same conclusion 
with certainty since there is a significant fraction of triads that may act unfavourably 
(see Merilees & Warn 1975). 



The preferential transfer of energy and enstrophy, together with the scale-selective dissi- 
pation by molecular viscosity, constitutes the backbone of the dual-cascade hypothesis. 
Let us consider the dynamical scenario described in the preceding subsection, with an 
additional assumption that the fluid is initially at rest. When the forcing region becomes 
unstable and the injected enstrophy gets transferred toward the high wavenumbers, palin- 
strophy P = (| VA-0| 2 )/2 is built up and growth of the enstrophy dissipation 2i/P ensues. 
Meanwhile growth of enstrophy is due solely to enstrophy injection and is suppressed by 
its own growing dissipation. If the growth of palinstrophy is sufficiently rapid, the en- 
strophy can quickly become steady, achieving a value sufficiently low that 2vZ -C e. This 
allows for a strong inverse energy cascade to be realized. By 'strong' we mean that the 
inverse cascade carries virtually all the energy injection to the large scales; otherwise 
the inverse cascade is 'weak'. At least a transient direct enstrophy cascade seems plau- 
sible since the palinstrophy is concentrated mainly at high wavenumbers, presumably at 
k > k v ^> s (the so-called dissipation range in the KLB theory). The realization of the 
classical k~ 3 enstrophy cascade as a persistent rather than just transient phenomenon 



3.1. The preferential transfer 




(3.1) 



(3.2) 



3.2. The dual- cascade hypothesis and Kraichnan conjecture 
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requires that the high concentration of palinstrophy around k v (which may be correctly 
termed a 'palinstrophy condensate') remain stable, withstanding the enormous dissipa- 
tion of palinstrophy in that region and not redistributing itself in wavenumber space. 
Note that if viscosity is replaced by a scale-independent dissipation, i.e. a mechanical 
friction, then both energy and enstrophy are dissipated at the same rate. The above 
argument for an inverse cascade fails to apply in this case. If the usual viscosity term 
vA 2 i() is replaced by an inverse viscosity represented by, for example, fitp, then the dy- 
namical behaviours of energy and enstrophy are reversed: energy behaves as enstrophy 
and vice versa. In this case the preferential transfer of energy and enstrophy gives rise to 
a predominant increase of the energy dissipation /Lt(|?/j| 2 ). The energy can then quickly 
become steady, achieving a value sufficiently low that 2/j.E <C s 2 e, i.e. the enstrophy 
dissipation is much less than its injection rate. This allows for a direct enstrophy cascade 
to be realized. Of course this case has no physical basis and is employed here only to 
illustrate the effects of the scale selectivity of the dissipation on the transient dynamics. 
We note in passing that numerical simulations aiming to verify the dual-cascade picture 
(or part of this picture: either an inverse or a direct cascade) employ both hyperviscosity 
and inverse viscosity. The results in this direction constitute a rich literature; some re- 
cent studies are|Paret, Jullien fc Tabeling (1999)[| Boffetta, Celani fc Vergassola ( 2000) 
Lindborg fc Alvelius (2000)| and |Chen et al. (2003)| The present arguments suggest that 



the effects of these scale-selective dissipation mechanisms on the transient dynamics (and 
probably beyond) should be taken into consideration when interpreting the results. 

For finite Reynolds numbers a small but non-negligible fraction of the injected energy 
is dissipated, resulting in a weak (i.e. less than complete) inverse energy cascade. For 
the fc -5 / 3 energy range this dissipation occurs mainly around the forcing wavenumber s, 
resulting in a similar fraction of enstrophy dissipation around s. Now if the remainder of 
the enstrophy injection were to be transferred to and dissipated at k 3> s, there would 
necessarily be a severe step in the spectrum between the forcing region and the rest of the 
enstrophy range. This means that for power-law scalings (without such a step), a weak 
direct enstrophy cascade is not permitted. This is in contrast to the robustness of the 
weak inverse cascade, which is readily observable in numerical simulations even in the 
complete absence of an accompanying direct enstrophy cascade (see Tran 2004 and Tran 
& Bowman 2004). Hence a plausible dynamical scenario is that for moderate Reynolds 
number the enstrophy dissipation occurs throughout the direct-transfer range up to some 
high wavenumber, which depends on the Reynolds number. The energy spectrum of this 
enstrophy dissipation range (instead of enstrophy inertial range) should then scale as 
fc~ 5 . The question is whether, for progressively higher Reynolds number, the inverse 
energy cascade can become stronger and a direct enstrophy cascade realizable. In a more 



quantitative analysis, Tran (2004)| suggests that a quasi-steady state featuring an inverse 
cascade of arbitrary strength (via the Kolmogorov-Kraichnan fc -5 / 3 spectrum) and a 
uniform dissipation of enstrophy among the wavenumber octaves of the direct-transfer 
range is plausible. This picture is consistent with the preferential transfer of energy 
and enstrophy, required by the conservation laws, and explains the numerical results 



of Tran (2004) and Tran & Bowman (2004) mentioned earlier. It may also explain the 



numerous numerical results targeted at the inverse energy cascade, for which there is 
hardly an enstrophy range due to limited resolution. 

In order to apply the dual-cascade hypothesis to turbulence confined to a bounded 



domain, Kraichnan (1967) suggests that after reaching the lowest available wavenumber 
fco, the inverse energy cascade maintains its strength during the quasi-steady stage, de- 
positing energy onto fco (or in the vicinity of fco, as presently considered). This process 
continues until the growth of energy at or around fco is limited by its own dissipation, 
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resulting in a huge pile-up of energy and enstrophy in this region of the spectrum. Note 
that even in the picture of no direct enstrophy cascade discussed in the preceding para- 
graph, the Kraichnan conjecture still applies although in this case the concentration of 
energy and enstrophy at the condensate would not be as dramatic as in the original 
dual-cascade case. This topic is discussed further in § 5. 

4. Inverse energy transfer 

In this section, we derive a rigorous upper bound for the energy that gets transferred 
across a low wavenumber I. It is shown that if a power-law scaling ak~ a is assumed for 
the energy spectrum in the inverse-transfer range and if a < 5/3, then no significant 
fraction of the energy injection can get transferred across £ for sufficiently low £. This 
result implies that the —5/3 slope represents the minimal steepness of the energy inertial 
range that can support an inverse energy cascade. For the special case I = kg a sharp 
estimate of the energy transfer onto ko is obtained. 

4.1. Inverse energy flux 

We assume that the spectral support of / is bounded from below by a wavenumber so 
and that the energy injection is bounded. The usual requirement of spectral localization 
of / can be relaxed for most of this section. For t < so the governing equation for the 
evolution of the large-scale energy density E < — (\Vip < \ 2 ) /2 is obtained by multiply- 
ing the governing equation 1|2.1[) by ip < and taking the spatial average of the resulting 
equation: 

^E K = (>>p < J(i>,AiP))-2 V Z < = -{^J{^,^<))-2vZ < , (4.1) 

where the second equality is due to (|2.2p . The nonlinear term represents the energy that 
gets transferred into the low- wavenumber region [k , £}, i.e. the energy flux across £, which 
drives the large-scale dynamics. 

We now derive an upper bound for the nonlinear term in l|4.1|) and then estimate this 
bound for power-law energy spectra. The steps go as follows 

|(AV>JW>, </><)) I = l^li,^)!^,^))! 

= \{i,J^>^<)+^>^<))\ 

< (i^iiv^liv^l + l^livviiv^l) 
<<|v^>||v^|(|v^<| 2 + |v<a<| 2 ) 1/2 ) 

<2^ /2 ^ sup( |vV,<| 2 + |V^<| 2 ) 1 / 2 

X 

= 2ElJ 2 E^ sup(|^<J 2 + |^< | 2 + 2|V< | 2 ) 1/2 , (4.2) 

where the first step is due to l|2.6|) and all the subsequent steps involve simple manip- 
ulations and the Cauchy-Schwarz inequality. Here E > = (|V'0 > | 2 )/2 is the 'small-scale' 
energy density associated with wavenumbers k > £. We observe that the upper bound for 
|AV' < | in (|2.9|) is also an upper bound for \ip^ x \, \4>y y \ and 2\ip< y \. Hence we can deduce 
that 

|(A^^,^)}| < c'f Z^E^E 1 ' 2 , (4.3) 
k 
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where d = 1\fic. The appearance of Z < in l|4.3[l is due mainly to l|2.6|l . which enables us 
to 'transfer' one spatial derivative from ip to if) < . In the velocity formulation of the NS 
system, the calculation leading to the upper bound l|4.3|) for the inverse flux is straight- 
forward. One of the advantages of the vorticity formulation, when equipped with (|2.6|) . 
is that various nonlinear terms, for example the term (A™?/; J(ip, A-0)) in the evolution 
equation of the quadratic quantity (ipA n+1 ip) /2, where n is an integer, can be manipu- 
lated with ease. Since Z < becomes smaller for progressively lower £, the right-hand side 
of (|4.3[) can remain bounded in the limit £ — * 0, even though both E > and E diverge in 
that limit. 

For further analysis of the right-hand side of l|4.3[) . we assume that an inverse-cascading 
range ak~ a has been established beyond I and possibly has reached ko, but no reflection 
or accumulation of energy has yet occurred. For this spectrum, E, E > and Z < can be 
estimated as follows: 



oc 



E < a k- a dk = hi' ' for a > 1, (4.4) 

Jko " - 1 

p CO 

E > <a / k- a dfc = — —I x ~ a for a > 1 (4.5) 



and 

Z < =al k} a dfc < — — — £^~ a for a< 3. (4.6) 

Jko 3 - a 

The requirement a < 3 is consistent with a persistent inverse energy cascade because a 
steeper spectrum would render the inverse-cascading range dissipative, which would not 
support such a cascade. The requirement a > 1 poses no loss of generality as becomes 
apparent shortly. By substituting the above estimates into (14.311 we obtain 



„/„3/2 / p \ (a + l)/2 

i(AM^))i< (a _;; 3 _ a)V2 (£) (4.7) 



Y - ] 

Now in the limit fcp — > 0, we have 



Z >a [ ° k 2 - a dk = -?—s?- a . (4.8) 
Jo 3 - a 

Solving for a from l|4.8|l and substituting the result into (|4.7|l we obtain 

KAM^^))! < \ 7T - ■ (4-9) 



(a-l)sg V^o/ V s o, 

An interesting conclusion can be readily drawn from l|4.H|l . Given a fixed ratio £/ko and 
finite enstrophy density Z (recall that we must have Z < e/2v), if a < 5/3, then the right- 
hand side of 1)4. 9|) can be made arbitrarily small provided £/sq is sufficiently small. In 
other words, no energy inertial range shallower than fc~ 5 / 3 would be capable of sustaining 
an inverse energy cascade that carries a nonzero fraction of the energy injection to the 
low-wavenumber region [fco,^], for sufficiently low t. Thus the Kolmogorov-Kraichnan 
spectrum represents the shallowest possible spectrum that is necessary for the 
existence of the classical inverse energy cascade. 



4.2. Energy transfer onto ko 

When £ = ko, the upper bound for the inverse flux derived in the preceding subsection 
becomes an upper bound for the energy that gets transferred onto fco- For this special case, 
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some improvement on the bound is possible and a sharp estimate for the Kolmogorov 
constant can be derived. 

We employ the notation ip — ip K + ip > as in the previous sections, where £ = ko- 
That means ip < , which is replaced by ip in what follows, consists of only four degenerate 
components corresponding to the four wavevectors (±fco,0) and (0, ±fco)- Namely, 

ip° — pi exp{jfcox} + ip\ exp{— ikox} + (p 2 exp{ikoy} + (p 2 exp{— ikoy}. (4-10) 

Inequality l|2.9ll becomes 

\AiP°\ < 2fc 2 (M + 1^1) < 2 3 / 2 fc (|fc Vi| 2 + k 2 \^\ 2 ) 1/2 = 2 3 / 2 k y 1 2 /2 (k a ), (4.11) 

where ^ 2 (ko) = (I Vip°\ 2 )/2 is the modal energy associated with ko. 

The energy that gets transferred onto ko, i.e. the nonlinear term (ip° J(ip, At/')), can 
be estimated as follows: 

|(^ J(^,A^))| = \(4,°JW>, AiP>))\ = \(Aip>J(ip>,ip°))\ 
= \(iP>J(iP>,^ x )) + (iP>J(iP>,^ y ))\ 

<2 3 / 2 ko* 1 2 /2 (k Q )E > , (4.12) 

where the replacement of ip by ip > = ip — ip° in the second step is a consequence of both 
Aip° = -k 2 ip° and (|2.2p. the third step is due to (|2.6I) . the last step is due to (|4.11[1 and 
all other steps are straightforward. Here E > = (|VV' > | 2 )/2 is the total energy density 
with the contribution from ko removed. The nonlinear term in the last equation of l|4.12|l 
represents an upper bound on the energy that gets transferred onto ko . 

As in the previous case we assume an inverse-transfer range ak~ a down to ko and 
estimate the energy that gets transferred onto ko, before the arrest of the inverse cascade 
would deform the assumed power-law spectrum. The factors of the nonlinear term in 
<|4.12[1 can be estimated as follows: 

r°° 2( 1 " Q )/ 2 o , 

E><a k~ a dk = kl~ a , 

JV2k " - 1 

* 2 (M = ^ = f (4.13) 

Z7T Z7T 

where 2 1//2 fco is the second lowest wavenumber, and E(ko) — ak^ a is the energy supported 
at ko- It follows that 

2^k Q ^\ko)E > < 2l ~ a/2a '[\ t 3a)/2 - (4-14) 
n{a — 1 j 

Like the argument in the preceding subsection, ineq. (|4.14|) together with l|4.12|l implies 
that an inverse-cascading range shallower than fc~ 5 / 3 is incapable of supporting a nonzero 
transfer of energy to ko for sufficiently low fco. 

4.3. Bounds for the Kolmogorov constant 

Suppose that the turbulence is driven by a steady energy injection rate e and that an 
inverse cascade carrying a fraction r of this injection toward fco y i a a afc -5 ^ 3 energy 
inertial range has been established. The left-hand side of 1|4.14(1 , which is an upper bound 
for the inverse energy transfer onto fco, cannot be smaller than re. Hence we can deduce 
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from (|4.14jl that 

a>{—) • (4.15) 

It follows that the Kolmogorov constant C which appears in the energy spectrum of the 
classical energy inertial range as E(k) = Ce 2 / 3 fc~ 5 / 3 is bounded from below by 

/2 5 / 6 7rr\ 2/3 

C>(^) - (4.16) 

Note that in the classical case r = 1. 

An upper bound for C can be derived on the basis of the classical spectrum alone. We 
observe that the energy dissipation by the energy range cannot exceed (1 — r)e. Hence 
we have 

2vCe 2/3 [ k 1/3 dk < (1- r)e. (4.17) 

It follows that 



o 



If we assume that e is injected around a forcing wavenumber s, so that the enstrophy 
injection is given by r\ = s 2 e, then we can rewrite l|4.18|l as 

C< { J V , (4.19) 

where r = 2^s 2 is the dissipation rate at the forcing wavenumber s. The upper bound 
H4.19fl is expected to be sharp in the limit of high Reynolds number. In fact l|4.19|) is 
essentially an equality in that limit since the energy dissipation by k > s is negligible. 
For a fixed enstrophy injection rate 77, constancy of C requires that the ratio (1 — t)/t 
remain constant. On the other hand, C may diverge in the limit of infinite Reynolds 
number (r — * 0) provided the inverse-cascade strength r approaches unity slower than r 
tends to zero. This dynamical scenario, which has been touched upon in § 3, is considered 



by Tran (2004) for the unbounded case, where it is suggested that an inverse energy 
cascade that is progressively stronger with progressively higher Reynolds number is quite 
plausible, requiring no boundedness of enstrophy or of any other quadratic quantities in 
the limit of infinite Reynolds number. In this picture there is no anomalous enstrophy 
dissipation: the enstrophy dissipation is uniformly distributed among the wavenumber 
octaves of the direct-transfer range, whose energy spectrum scales as fc -5 . 



5. Power-law spectra vs. Kraichnan's condensate 

In this section we examine the Kraichnan picture for steady dynamics in a bounded 
domain, and derive constraints on the enstrophy and palinstrophy distribution for such a 
picture. We then contrast this picture with power-law spectra. For simplicity we follow the 
non-average treatment in the preceding sections. Nevertheless all dynamical quantities, 
including the energy spectrum, can be understood in an appropriate average sense since 
such a reinterpretation requires a straightforward reformulation of the problem. 

By virtue of the conservation laws, the evolution of the energy and enstrophy is gov- 
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erned by 

dE 

— = -2vZ + e, (5.1) 

< ¥i = -2vP + r 1 , (5.2) 

where e = —{ftp} and rj = (fAtp) are, respectively, the energy and enstrophy injection 
rates. The usual assumption of spectral localization of / around s is invoked, i.e. rj = s 2 e. 
The forced-dissipative balance of both energy and enstrophy implies (see Constantin et 
al. 1994, Eyink 1996, Tran & Shepherd 2002 and Kuksin 2004) 

P = s 2 Z, (5.3) 

or in terms of the energy spectrum E(k), 

(fc 2 - s 2 )k 2 E{k) dk = 0. (5.4) 

Equation (|5.4|) is the focus of the present section. 



5.1. Energy and enstrophy condensate 

Eq. 1)5. 3|1 . which is independent of Reynolds number, has been interpreted as implying 
that the enstrophy dissipation must occur in the vicinity of s, thus precluding a direct 
cntrophy cascade. However, this interpretation assumes that all the enstrophy is actively 
involved in the enstrophy cascade. As noted earlier, the Kraichnan energy condensate 
is also an enstrophy condensate, in which case most of the enstrophy is trapped in the 
condensate and only a tiny amount is free to participate in the enstrophy cascade. Thus, 
this interpretation implicitly excludes the Kraichnan scenario. 

In fact, for spectra such that Z is distributed mainly around kg and P is distributed 
mainly at k ^> s, the balance l|5.3|) can still hold (see below). Another example is that P 
has an equal contribution from each of the wavenumber octaves [s, 10s], [10s, 10 2 s], • • •, 
up to some wavenumber k v « s 2 /fco, provided that Z has a similar contribution from 
the wavenumber octaves of [ko, s]. This uniform distribution of the energy and enstrophy 
dissipation in their respective ranges (requiring k energy and A; -5 enstrophy ranges) 



is considered by Tran & Bowman (2003) and Tran (2004) (also see below). The former 
example includes the Kraichnan scenario, for which the energy of the condensate E < (with 
£ « ko) satisfies E< i=s e/2vk^. The energy dissipation occurs mainly at the condensate 
since 2vZ < w e, where Z < fc§i?< is the enstrophy of the condensate. The enstrophy 
dissipation by the condensate is 2vP < w k 2 e, where P< « k 2 Z< is the palinstrophy of 
the condensate. This dissipation is negligible as compared with the enstrophy injection 
s 2 e if k /s <§; 1, i.e. if the energy inertial range is sufficiently wide. This possibility allows 
for the spectral distribution of palinstrophy, i.e. the spectral distribution of enstrophy 
dissipation, to be concentrated at the high wavenumbers, while posing no threat to the 
balance P = s 2 Z. 

The discussion in the preceding paragraph implies that the Kraichnan energy conden- 
sate would also be an enstrophy condensate. The extent of this highly localized enstrophy 
distribution at the condensate would be comparable to that of the palinstrophy around 
the dissipation wavenumber k v . To see this explicitly let us rewrite the balance equation 
l|5.3|) as 

P<+P> =s 2 (Z < +Z > ), (5.5) 



where the subscripts '<' and '>' refer to the 'trapped' (in the condensate) and 'free' 



Remarks on the KLB theory of 2D turbulence 



11 



quantities. It follows that 

where k 2 — P < /Z < w kg. In the Kraichnan picture the ratio P > /s 2 Z > ss k 2 j 1 s 2 diverges 
in the limit of infinite Reynolds number. Hence the ratio Z < /Z > diverges in a similar 
manner. Thus the Kraichnan energy condensate is also an enstrophy condensate. 

5.2. Power-law spectra 

Following Constantin et al. (1994), Tran (2004) and Tran & Bowman (2003) we assume 
a two-range energy spectrum: 



where a and (3 are constant and k v marks the end of the fc - ' 3 range, beyond which the 
spectrum becomes steeper. This approximation a priori excludes the Kraichnan conden- 
sate, but allows for some 'milder' accumulation of energy at ko. For example, the case 
a = 3 previously considered would allow for most of the system's energy to reside in 
a few lowest wavenumbers. Substituting l|5.7(l into l|5.4J) allows one to derive interest- 
ing constraints on a and (3. By rewriting (|5.4J) in the form f^(s 2 — k 2 )k 2 E(k) dk = 
J^°(k 2 — s 2 )k 2 E(k) dk, using the approximation l|5.7|l for E(k) and making the substi- 
tutions k = k/s for k < s, and k — s/k for k > s, |Tran &: Bowman (200*3) show that 



(1- n 2 )K, 2 - a dK,> / (1 - k 2 )k,P- 6 d/s, (5.8) 

ko / s J sj k v 

where the inequality is a consequence of dropping from the equality the contribution 
beyond k v , i.e. the quantity D v = J^°(k 2 — s 2 )k 2 E(k) dk. Except for the approximation 
H5.7J1 . the constraint (|5.8|l is rigorous. In the spirit of KLB {k v — » oo as v — > 0), we 
consider the case kg/s > s/k v , which corresponds to high- Reynolds number turbulence. 
It is then easy to see from (|5.8[) that 2 — a < (3 — 6, or equivalently, 

a + (3 > 8. (5.9) 

For moderate Reynolds number f3 > 5 (see Tran 2004 and Tran & Bowman 2004), 
so we first consider the case (3 > 5. The quantity D v is then negligible, making l|5.8|l 
essentially an equality and l|5.9|l essentially an equality if k^/s ~ s/k v . In this case we 
obtain (a, (3) — (3 — 6, 5 + S), where the limit 6 — > is expected for high Reynolds 
numbers. This solution with 8 = means that the energy and enstrophy are uniformly 
dissipated among the wavenumber octaves of their respective ranges. On the other hand, 
if j3 < 5 (which cannot be logically ruled out, just as the Kraichnan picture cannot be a 
priori excluded), then a > 3. In this case the inequality ">" in (|5.8|) could become "^>" 
since D v could become large. As a consequence a slight decrease of (3 from the critical 
value (3 = 5 must be met with a much more significant increase of a from a = 3. This 
means that a could quickly approach (3 before the latter would drop significantly below 
its critical value. Hence the two-range steady spectrum Ij5.7ll . when supplemented by the 
usual condition a < (3, obeys the constraint (3 > 4 + S, where 6 > depends on how 
E[k) decays for k > k v . The case (3 < 4 is inadmissible and would necessarily require 
an energy condensate. For the Kolmogorov-Kraichnan energy inertial range (a — 5/3), 
a Kraichnan-type condensate would also be required even when [3 = 5 since the integral 
on the left-hand side of (|5.8(l is less than unity (for all ratios ko/s) while its counterpart 



12 



C. V. Tran and T. G. Shepherd 



on the right-hand side is « ln(fc„/s), which diverges as k v /s — ► oo. Finally, the classical 
fc -3 direct-cascade range requires a k~ 5 inverse-cascade range if only power-law spectra 
are permitted. The enstrophy content of the Kraichnan condensate would be equivalent 
to the enstrophy content of a k~ 5 energy range. 



6. Conclusion and discussion 

In this paper the inverse energy transfer in 2D NS turbulence is studied. The main 
result obtained is an upper bound for the inverse energy flux. This upper bound is then 
estimated using power-law spectra. It is shown that a steady inverse energy flux requires 
that the energy spectrum of the inverse-transfer range be at least as steep as fc -5 / 3 : 
a shallower spectrum is unable to support an inverse energy cascade, and the inverse 
energy flux necessarily diminishes as it proceeds toward sufficiently low wavenumbers. 
This result provides for the first time a rigorous basis for the Kolmogorov-Kraichnan 
fc -5 / 3 energy inertial range in 2D turbulence. Other results include estimates of the 
Kolmogorov constant, a relation between the enstrophy and palinstrophy distribution 
for the Kraichnan condensate, and an analysis of power-law scaling in contrast to the 
Kraichnan conjecture. 

The inverse energy flux across a wavenumber £ in the energy range is found to depend 
on the system's energy and on the enstrophy content of the wavenumber region lower 
than I. The former is supposed to grow without bound in the limit of unbounded do- 
main, while the latter becomes progressively smaller for progressively lower I. These two 
effects counterbalance each other in such a way that the classical inverse energy cascade 
necessarily requires that the energy spectrum of the inverse-transfer range be at least as 
steep as the Kolmogorov-Kraichnan fc -5 / 3 spectrum. 

The derived upper bound for the inverse energy flux has been used to deduce a 
lower bound for the Kolmogorov constant C in the classical energy spectrum E{k) = 
Ce 2//3 fc -5 / 3 . This bound is of order unity and is expected to hold whether or not C is a 
universal constant. 

We have elaborated on the Kraichnan conjecture of energy condensation. A constraint 
on the enstrophy and palinstrophy distribution for this picture has been derived and 
contrasted with those for power-law scalings. For power-law spectra, a fc -3 energy range 
and a k~ 5 enstrophy range are consistent with the forced-dissipativc balance of energy 
and enstrophy. This new picture has some justification from the numerical results of 
Borue (1994) Tran (2004) and Tran fc Bo wman (2004) at least for moderate Reynolds 
numbers. A Kraichnan-type condensate is required for energy-range spectra shallower 
than fc~ 3 even if the enstrophy-range spectrum remains as steep as fc~ 5 . The energy and 
enstrophy level of the condensate would dramatically increase should both the energy 
range be shallower than fc -3 and the enstrophy range be shallower than fc -5 . The classical 
fc -3 enstrophy range would require an excessively high concentration of energy and en- 
strophy at the condensate: the enstrophy content of the condensate would be equivalent 
to that of a fc -5 energy range. 

Several authors have considered steady bounded turbulence and, on the basis of the 
balance equation H5.3(l . ruled out the existence of a direct enstrophy cascade (see Con- 
stantin et al. 1994; Eyink 1996, p. 110; Tran & Shepherd 2002; Kuksin 2004). These 
works either explicitly or implicitly exclude a priori the Kraichnan conjecture in favour 
of power-law spectra. In particular Constantin et al. (1994) assume the power-law scaling 
described by l|5.7(l with the classical exponents a = 5/3 and (3 = 3, and find that this 
spectrum is incompatible with the balance P = s 2 Z. However, the two-range spectrum 
of KLB is supposed to be only quasisteady, with the inverse energy cascade proceeding 



Remarks on the KLB theory of 2D turbulence 



13 



toward ever-lower wavenumbers, carrying with it virtually all the energy injection. In 
other words, the KLB spectrum is supposed to correspond to dE/dt = —2vZ + e e, 
or equivalently, to P > s 2 Z, not to the balance P = s 2 Z, which is achieved only in 
a steady state. The quasisteady limit as t — > oo is clearly incompatible with a finite 
domain. Eyink (1996)| and |Tran fc Shepherd (2002) infer from 15.3(1 that the dissipation 
of enstrophy occurs mainly in the vicinity of the wavenumber given by ^/P/Z, which 
is just the forcing wavenumber s for steady turbulence, thereby ruling out the existence 
of a direct enstrophy cascade. This conclusion is valid for moderate ratios s/kg, i.e. for 
relatively narrow energy ranges, and for various 'regular' spectra such as ones with the 
energy range shallower than fc -3 . However if the Kraichnan conjecture holds, then this 
conclusion would not be valid. Kuksin (2004) employs a special forcing, such that both 
the energy and enstrophy injections are proportional to the viscosity coefficient v, to 
show that in the limit of infinite Reynolds number the statistical equilibrium palinstro- 
phy remains bounded. This result leads Kuksin to conclude that the KLB theory does not 
apply to the statistically steady dynamics of bounded NS turbulence in a doubly periodic 
domain. Similar to |Eyink (1996)| and |Tran fc Shepherd (2002)| the balance P = s 2 Z is 
also obtained, but no further constraints related to the Kraichnan scenario are derived. 
The fact that P remains bounded in the limit of infinite Reynolds number is due to the 



special choice of the forcing. In essence this result is equivalent to that of Eyink (1996) 



and Tran & Shepherd (2002) except by a scaling factor of the inverse of the Reynolds 
number. Hence the above analyses apply to this case as well. In particular the Kraich- 
nan conjecture of energy condensation cannot be logically excluded. The present work 
emphasizes the need for some treatment beyond the forced-dissipative balances of the 
energy and enstrophy, before such an exclusion can be fully justified. 
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